Draft version February 5, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 6/22/04 



REDSHIFTS FROM SPITZER SPECTRA FOR OPTICALLY FAINT, RADIO SELECTED INFRARED 

SOURCES 

D. W. Weedman 1 , E. Le Floc'h 2 , S. J. U. Higdon 1 , J. L. Higdon 1 , J. R. Houck 1 

Draft version February 5, 2008 

ABSTRACT 

Spectra have been obtained with the Infrared Spectrograph on the Spitzer Space Telescope for 18 
optically faint sources (R > 23.9 mag) having f„ (24/im) > 1.0 mJy and having radio detections at 20 
cm to a limit of 115 fj,Jy. Sources are within the Spitzer First Look Survey. Redshifts are determined 
for 14 sources from strong silicate absorption features (12 sources) or strong PAH emission features 
(2 sources), with median redshift of 2.1. Results confirm that optically faint sources of ~1 mJy at 
24 /xm are typically at redshifts z ~ 2, verifying the high efficiency in selecting high redshift sources 
based on extreme infrared to optical flux ratio, and indicate that 24 /jm sources which also have radio 
counterparts are not systematically different than samples chosen only by their infrared to optical 
flux ratios. Using the parameter q = log[fi,(24 /im)/f„(20 cm)], 17 of the 18 sources observed have 
values of 0<q<l, in the range expected for starburst-powered sources, but only a few of these show 
strong PAH emission as expected from starbursts, with the remainder showing absorbed or power-law 
spectra consistent with an AGN luminosity source. This confirms previous indications that optically 
faint Spitzer sources with f„ (24/zm) > 1.0 mJy are predominately AGN and represent the upper end of 
the luminosity function of dusty sources at z ~ 2. Based on the characteristics of the sources observed 
so far, we predict that the nature of sources selected at 24/im will change for f l/ (24/im) < 0.5 mJy to 
sources dominated primarily by starbursts. 

Subject headings: dust, extinction — galaxies: high-redshift - radio: galaxies X-ray: galaxies — 

infrared: galaxies — galaxies: starburst — galaxies: AGN — 



1. INTRODUCTION 

Mid-infrared imaging surveys underway with the 
Spitzer Space Telescope demonstrate the presence of 
large numbers of sources whose luminosity must arise pri- 
marily from dust re-emission and which are significantly 
obscu red at optical wavelengths (e.g. iPapovich et al. I 
( 2004)). The source counts are consistent with expec- 
tations derived from efforts to account for the total in- 
frared background and have been modeled as showing 
the evolution of lum in ous, star-forming g alaxies (e.g. 
lLagache et a!~l l|2004|) . iCharv et all l|2004in . Spectro- 
scopic and photometric redshifts also indicate that the 
majority of such sources have in frared luminosity pow- 
ered primarily by star formation l|Perez-Gonzalez et al. I 
2005). Confirming the nature of these sources and 
their redshift distribution is crucial to understanding 
the evolution of star formation and AGN activity in 
the universe, especially at early epochs. A similar sit- 
uation has existed for much longer in efforts to under- 
stand the nature of o p tically faint radio sourc es (e.g. 
iRichards et aT~l (fl99flh . iFomalont et~aT~l pOOl ). The 
observation that many sub-mJy radio sources are iden- 
tified with a population of faint, blue galaxies initially 
indicated that most optically-fain t radio sources are pow - 
ered primarily by star formation l)Haarsma et al. 1 (2000') . 
This had also been indicated by the correlation be- 
tween infrared and radio fluxes for starb u rst systems 
llCondon et al. I (|1982fl . iHelou et al. I l|1985lUYun et al. I 
(2001)). However, recent studies have shown that many 
faint (~ 100 /iJy) radio sources do not show the infrared 
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fluxes llH igdon et al. ll20 05) or submillimeter detections 
l|Chapman et al. 1120041) expected for starburst systems. 
It is important to determine, therefore, what fraction of 
the faint radio source population does have firm indica- 
tions of starbursts. 

A major observational challenge is to obtain redshifts 
or spectral diagnostics for infrared and radio sources 
too faint for optical spectroscopy, sources having opti- 
cal magnitudes > 24 mag. For sources that are suffi- 
ciently bright in the infrared (f^ (24/im) > 0.75 mJy), 
redshifts can be determined to z ~ 2.8 using the Infrared 
Spectrograph on Spitzer (IRS), based on strong spectral 
features from silicate absorption or PAH emission. Re- 
sults for the first set of sources selected only on the ba- 
sis of extreme IR/optical flux ratios indicated that such 
sources were typically at z ~ 2 and were usually simi- 
lar to the absorption spect ra of local sources powered by 
AGN (jHouck et all l|2005D : hereinafter H05). However, 
some s ources selected u sing colors characterising star- 
bursts l| Y an et a f~l I2005T) or using submillimeter detec- 
tions (|Lutz et al. I2005|) showed PAH emission character- 
istic of starbursts. Submillimeter observations of Spitzer 
24/im sources selected because of extreme IR/optical ra- 
tios did not show the submilli meter fluxes expecte d for 
a star burst-dominated sample ijLutz et al. Il2005blk in- 
stead, most of these sources typically have stronger mid- 
infrared fluxes, consistent with the presence of hotter 
dust powered by an AGN. It is clear, therefore, that 
the optically faint, infrared selected population has both 
AGN and starburst constituents, and further samples 
are needed to define their redshift distributions and rela- 
tive fractions in the overall luminosity function of dusty 
sources. 

The Spitzer First Look Survey has imaged 4.4 deg 2 at 
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24 fim with the Multiband Imag ing Photometer (MIPS) 
instrument llRiek e et al. | 2004fl: initial results are de- 
scribed by IMarleau et al. I l)2004|) although no catalogs 
are as yet available. The survey also has accompany- 
ing deep imaging with the Very L arge Array (VLA) for 
which a source catalog is available {Condon et al. ll2003ft 
and imaging in R band from the National Optical As- 
tronomy Observatory 3 . The initial comparison of 24 /im 
flux densities with 20 cm flux densities (a ratio parame- 
terized by q = log[f„(24 /xm)/f !/ (20 cm)]) indicated that 
sources which are detected in both in frared and radio 
cluste r around a median q value of 0.8 {Appleton et al. I 
2004), which is the ratio expected from the previously 
known radio- infrared correlations for starbursts. This 
apparent agreement with expectations is misleading, be- 
cause it does not include the large numbers of sources 
detected in the radio but not in the infrared, or vice- 
versa. A comparison of Spitzer and VLA 20 cm surveys 
to similar detection limits, although in a much smaller 
sky area, indicated that only 9% of the 24 /im sources are 
detected in the radio, and only 33% of the rad io sources 
detected in the infrared {Higdon et al. 1 12005^> . Taking 
into account these limits, radio sources have a median q 
that is negative whereas infrared sources have a median 
q of a bout unity. This result was used bv lHigdon et al. 1 
( 2005) to conclude that the majority of radio sources in 
these VLA surveys are powered by AGN, whereas the 
majority of infrared sources are powered by starbursts. 

The result of starburst dominance for the Spitzer 
24 /im sources based on various survey comparisons and 
modeling of counts is inconsistent with the conclusion of 
H05 based on IRS spectroscopy that most optically faint 
sources of ~ lmJy at 24 /im are distant ultraluminous 
galaxies powered by AGN. They observed 31 sources in 
the Boote s field of the NOAO Dee p Wide-Field Survey 
fNDWFS. ITannuzi and Devi {199^ having R > 25 mag 
and i v (24/im) > 0.75 mJy. Of the 17 sources with deter- 
minable redshifts, 16 were best fit with heavily absorbed 
templates and have median z of 2.2. Because the ab- 
sorbed spectra for the local templates arise from objects 
with AGN (the prototype being Markarian 231), and 
because of the radio properties and luminosities of the 
sources, H05 interpret the heavily absorbed sources as 
bring primarily powered by obscured AGN. S pectroscopy 
of 8 s ources in the Spitzer FLS selected bv lYan et al. I 
(2005) using color criteria targeted to select starbursts 
showed that 2 of the 6 sources with measurable spec- 
tral features show s trong PAH e missio n features. Two 
sources observed bv lLutz et al. 1 {2005|) selected because 
of previous SCUBA detections also show PAH emission. 
It is clear, therefore, that the optically faint infrared pop- 
ulation contains a variety of sources, and selection crite- 
ria other than simply the infrared to optical ratio need to 
be invoked in efforts to sort and classify this population. 

To extend our samples and to begin investigation of 
the nature of infrared/radio sources, we undertook spec- 
troscopic observations of sources in the FLS area selected 
with infrared and optical criteria similar to those in our 
Bootes survey but with the additional criterion of a radio 
detection at 20 cm. 

3 FLS data products are available at 
http: / / ssc. spitzer . caltech.edu /fls 



2. OBSERVATIONS AND DATA ANALYSIS 

Our selection criteria for spectroscopic targets were to 
use the MIPS and VLA surveys with the R band optical 
image to inspect all sources in the V LA 20 cm catalog 
of the FLS area {Condon et al. Il2003|) having i v (24/im) 
> 0.75 mJy. To do this, we obtained the 24 /im images 
from the Spitzer archive and produced an FLS catalog 
of 24/im sources, reduced wit h the MIPS Data Analy- 
sis Tool {Gordon et al. ll2005|) . Point source extraction 
was performed using an empirical point spread function 
(PSF) constructed from the brightest objects found in 
the 24/im image and fitted to all the sources detected 
in the data. The flux of each object is derived using 
the scaled fitted PSF after applying a slight correction 
to account for the finite size of the modeled point spread 
function. For the present study, we are only interested in 
examining sources with f„ (24/im) > 0.75 mJy whereas 
the FLS 5 sigma limit is about 0.15mJy; our sample is 
complete, therefore, for these bright sources. Our 24/im 
catalog was compared to the VLA 20 cm catalog, and all 
sources having VLA detections (brighter than 115/xJy) 
were selected for comparison to the NOAO R band cata- 
log, also available in the Spitzer archive, and all of the ra- 
dio/infrared sources having R > 23.9 mag or being opti- 
cally unidentifed were examined. (Our criteria for source 
association were that cataloged source coordinates agree 
to within 2 " .) Our final sample consists of the 26 sources 
in the FLS having f„ (24/tm) > 1.0 mJy which are in the 
20 cm catalog and have R > 23.9 mag. We have ob- 
tained spectra for 18 of these sources with the IRS; the 
remaining 8 sources meeting the se criteria were alr eady 
within the program described bv lYan et al. I {2005D and 
not accessible to us. 

Observations and results for the sources discussed in 
this paper are summarized in Table 1. Coordinates 
listed are the 24 /im coordinates, which were used for 
the IRS targets. Spectroscopic observations were made 
with the IRS Short Low module in order 1 only (SL1) 
and with the Long Lo w module in orders 1 and 2 (LL1 
and LL2), described in lHouck et al. 1 1)2004() 4 . These give 
low resolution spectral coverage from ^8 /im to ~35 /im. 
Sources were normally placed on the slits by offsetting 
from nearby 2MASS stars; in a few cases with no suf- 
ficiently nearby 2MASS stars, direct pointing without 
offsets was used successfully. 

Sources which are observed are only a few percent of 
the brightness of the background flux, so the background 
is the dominant source of noise, and accurate background 
subtraction is essential before extracting a source spec- 
trum. Background subtraction for LL1 and LL2 modules 
can utilize the background observed in the off-source or- 
der; e.g., the LL1 slit provides an observation only of 
background when the source is in the LL2 slit, and this 
background can be subtracted from the observation when 
the source is in the LL1 slit. All images when the source 
was in one of the two nod positions on each slit were 
coadded to obtain the image containing the source spec- 
trum for that nod position. The background which was 
subtracted from this coadded source-spectrum image in- 
cluded coadded background images of both nod positions 

4 The IRS was a collaborative venture between Cornell Univer- 
sity and Ball Aerospace Corporation funded by NASA through the 
Jet Propulsion Laboratory and the Ames Research Center 
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when the source was in the other slit and including im- 
ages from the other nod position in the same slit. This 
means that the background observation subtracted from 
a source observation includes three times the integration 
time as for the source. This improves signal to noise 
in the subtracted background. We experimented with 
the addition of background images from observations of 
different sources in an attempt to produce a standard 
background ("supersky") for application to all sources, 
but the background observed from source to source was 
not sufficiently stable to produce a reliable supersky. The 
differenced source minus background image for each nod 
position was used for the spectral extraction, giving two 
independent extractions of the spectrum for each LL or- 
der. These two were compared to reject any highly outly- 
ing pixels in either spectrum, and a final mean spectrum 
was produced. For SL1, there was no separate back- 
ground observation with the source in the SL2 slit, so 
background subtraction was done by differencing coad- 
ded images of the two nod positions in SL1. 

Data were processed with version 11.0 of the SSC 
pipeline, and the bed. fits files were used for our spec- 
tral extractions. Extraction of sou rce spectra was don e 
with the SMART analysis package ijHigdon et al. |[2 0041; 
location of spectra on the slit for extraction was done by 
individual examination because in several cases serendip- 
itous sources also fell on the slit, and adopted back- 
grounds had to exclude these sources. Typical extraction 
widths shown in the cross-dispersion profile displayed by 
SMART were 3 to 3.5 pixels, to minimize noise from pix- 
els containing little source flux. While use of an extrac- 
tion window this narrow reduces the source flux some- 
what compared to the calibration observations that uti- 
lize a wider extraction, the 24 fim MIPS fluxes are avail- 
able for determining normalization of extracted spectral 
fluxes. It can be seen from comparison of the spectra 
shown in the Figures with the MIPS fluxes in Table 1 
that the extracted spectra typically agree at 24 /im to 
within 10 % of the MIPS flux. Final spectra for analy- 
sis of redshifts and features were smoothed to approxi- 
mately a resolution element, applying boxcar smoothing 
of 0.2 fim for SL1, 0.3 /im for LL2, and OA fim for LL1. 
These resulting spectra are shown in Figures 1, 2 and 3. 

3. DISCUSSION 

Even in spectra of faint sources with poor S/N, spec- 
tral features with sufficiently large equivalent widths ex- 
ist in the mid-infrared spectra for redshift determina- 
tion. At one extreme, these features are the strong 
PAH emission features characteristic of starburst galax- 
ies. The other extreme shows strong absorption fea- 
tures, the strongest being silicate absorption, with no 
PAH emission. Examples of PA H spectra for galax- 
ies observed with the IRS are in iBrandl et al. (2004) 
for NGC 7714, in IBrandl et al~l <l20Q5j) for the mean 
of starbursts, and in iWeedman et al. I l)2005j) for NGC 
3079. The prototype absorbed AGN is Markarian 231, 
which shows an apparent excess above the continuum 
centered at 8 /im (rest frame) , caused by absorption on 
either side, followed by a deep trough of silicate absorp- 
ti qn centered at 9.7um; t he IRS spectrum is discussed 
in IWeedman et aT~l <|2005f) . An even more extremely ab- 
sorbed source is IR AS F00183-7111, with IRS spectra in 
iSpoon et al.l ((2004) . An intermediate template with both 



absorption and emission is Arp 220, for which we utilize 
unpublished IRS spectra. This is characterized primarily 
by absorption, but the excess at 8 /jm is accentuated by 
the presence of PAH 7.7/zm emission. 

Redshifts can be determined by seeking either the 8 /im 
"excess" and following absorption, or the set of strong 
PAH emission features. The strongest PAH feature is 
at 7.7 fim (rest frame), so a similar redshift would be 
derived even if a spectrum is ambiguous as to whether 
the strongest feature is the 8 fim apparent excess or true 
PAH emission; the interpretation of the source, however, 
would be very different with the two alternatives. PAH 
features of about half the strength of 7.7 fim are at 6.2 /im 
and 11.2 /im. In order to interpret a feature as PAH, 
we require an indication that at least one of the other 
PAH features is present with the correct shape and rel- 
ative flux as scaled to the 7.7 /im feature, although the 
11.2 /tm feature is often redshifted out of our observed 
spectral range. Redshifts are estimated in two ways: in 
the first, we simply seek the features described above 
and assign a redshift based on the peaks in the spec- 
trum. In the second, we fit a selection of templates that 
attempt a formal chi-squared fit to the full spectrum so 
that strengths of the features relative to the continuum 
are also considered. These templates are a pure PAH 
emission spectrum such as M82 or NGC 3079, Arp 220 
(with both PAH emission and deep silicate absorption), 
Markarian 231, and F00183-7111. The redshifts derived 
with both techniques are listed in Table 1. Based on 
the differences in derived redshift depending on the tem- 
plate used, and on the differences in redshift measures 
from fitting templates and from seeking only the indi- 
vidual features, we estimate that the uncertainty in red- 
shift is typically 0.1. The final classification of an object 
as characterized by PAH emission (em) or silicate ab- 
sorption (abs) is adopted as listed in Table 1. Spectra 
with a straightforward choice of best-fitting template and 
spectra for which no features could be found are shown 
in Figure 1. Spectra having definite features but with a 
more uncertain assignment of template fits are in Figures 
2 and 3. 

Because we are attempting to use the presence of 
PAH emission or silicate absorption as an indicator of 
starburst-derived luminosity or AGN-derived luminosity, 
respectively, it is important to determine which charac- 
terizes the spectrum even if the redshift is not affected 
by this choice. This is not always an unambiguous de- 
cision, and we present in Figures 2 and 3 examples of 
the 7 sources in which either interpretation (all emis- 
sion and no absorption, or all absorption and no emis- 
sion) can provide a possible fit to the spectral features 
seen. We present these examples primarily as illustra- 
tions of the cautions that must be used in interpreting 
these IRS spectra with poor S/N; it is important not 
to assign high confidence to features that may not be 
real. Some of this ambiguity may arise because sources 
actually are composite and would formally best be fit 
with varied combinations of PAH emission and silicate 
absorption spectrum. We do not feel, however, that the 
S/N is adequate for this to improve our results regard- 
ing the classification of absorption dominance or emis- 
sion dominance, but believe instead that the ambiguities 
arise primarily because of the difficulty in judging if a 
particular spectral peak is a real feature or is noise. The 
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ambiguous sources which we illustrate are numbers 1, 
5, 9, 11, 14, 17, and 18; their spectra and the alterna- 
tive fits are shown in Figures 2 and 3. The choice of 
fit (absorption or emission) for these objects depends on 
deciding whether seemingly strong spectral features are 
real rather than accidents of noise. We display all spectra 
to the longest wavelength, 35 /im, at which some signal 
may be meaningfully recorded for these faint sources. We 
caution, however, that any large features beyond 33 /im 
should receive very little weight, because noise spikes are 
often found at these longer wavelengths where the de- 
tector sensitivity is falling rapidly. As examples of the 
uncertainties introduced in the fits, sources 1 and 5 in 
Figure 2 could be well fit by PAH features if the spike 
at ~ 33 /im is a real feature, but if it is not real, the fit 
cannot be PAH. Conversely, source 11 in Figure 2 can- 
not be PAH if the collection of spikes between 29 /im and 
33 /im represents real signal. Our final decision on which 
fit to accept from these ambiguous cases is based on our 
overall judgment of which spectral features are real, but 
we illustrate the alternative fits so that the reader may 
make their own judgment and use these examples for 
comparison to future ambiguous spectra. There is only 
one case where the ambiguity in fit would result in a sig- 
nificant difference in redshift, because of ambiguity re- 
garding identification of features. This is for source 17 in 
Figure 3. This is the poorest fit we have, and assigning 
a redshift based on the Markarian 231 template requires 
ignoring as noise an apparent excess of flux near 31 /im. 
If this is considered as a real feature, it is consistent with 
a 7.7 /im feature at a higher redshift than any source we 
have yet found, and which poorly explains the continuum 
flux at shorter wavelengths. We illustrate source 17 as 
the example of the one source for which we cannot decide 
if it has a measurable redshift or not. In our final clas- 
sifications summarized in Table 1, we do not include a 
redshift for source 17 and provide only the best fit power 
law for this source. The fits we finally adopt for the re- 
maining ambiguous sources are given in Table 1; these 
fits are judged by whether the pure absorption or pure 
emission templates shown in Figures 2 and 3 provide the 
better overall fit to the observed spectra. Only two of 
the 7 ambiguous cases in Figures 2 and 3 are judged to 
be better fit by PAH emission whereas the remainder are 
assigned to the absorption fit, except for number 17 with 
an indeterminate redshift. 

For the total set of 18 sources, we are able to assign 
redshifts to 14. The median redshift of our sample is 2.1 
and that of H05 is 2.2. The most notable difference in 
the results for the present sample compared to H05 is 
that we find proportionally a few more sources which are 
dominated by PAH emission spectra. Two of 14 sources 
(numbers 3 and 9) are well fit with a pure PAH spec- 
trum, whereas only one of the 17 sources in H05 was fit 
by the PAH-dominated spectrum of NGC 7714. We have 
3 of 14 sources best fit by Arp 220, which contains PAH 
emission in the template, and there were also 3 of these 
in H04. This slight excess of PAH sources is the only 
difference between the present sample and that of H05 
which might be attributable to the radio selection, pre- 
sumed to favor starbursts. The majority of sources with 
redshifts (12 of 14) are characterized by templates dom- 
inated by absorption (including the Arp 220 template), 
as were 16 of the 17 in H05. For further discussion and as 



a working hypothesis, we assume a simple classification 
into AGN-powered sources or starburst-powered sources 
determined only by whether the source shows strong ab- 
sorption or strong PAH emission, respectively. There are 
4 sources without detectable features; for these, the index 
a of the power law which best fits the observed spectrum 
is given in Table 1. The classification we assume for these 
objects is that they are also AGN characterised by hot 
dust but at unknown redshift; the absorption may be too 
weak to detect, or may be at z > 2.8, where it would be 
out of the detectable range of the IRS spectra. 

It is important to note an important caveat regard- 
ing the classification of AGN sources. This arises be- 
cause of the metal-poor, compact starburst SBS 0335- 
052, This has proven to have a starburst spectrum in the 
mid-infrared which is unique among all of the starbursts 
observed with the IRS. Unlike all other galaxies observed 
because they were previously classified as starbursts 
llBrandl et al. II200E ). it has no indication of PAH fea- 
tures l)Thuan et al. I l|1999ftlHouck et al. I <|2004bt^ and 
is the most extreme in this respect of all bl ue compact 
dwar f starbursts so far observed with the IRS l|Wu et al. I 
2005). It also shows weak silicate absorption. As a result, 
its observed mid-infrared spectral characteristics are very 
similar to the characteristics which we assign to AGN- 
powered sources. Despite its unusua l spectral character - 
istic s compared to other starbu rsts. iHunt et al. I (2005) 
and iHunt and Maiolinol l)2005|) suggest that SBS 0335- 
052 is the best local example for the spectral shape of 
a primordial starburst. If this is so, and if such ob- 
jects exist at luminosities and redshifts comparable to 
the sources we have found, the fraction of sources which 
we classify as AGN-powered is smaller than we have es- 
timated. 

Adopting the AGN classifications as described, it 
means that our sample has either 16 of 18 AGN, if 
Arp 220 absorption templates are assigned as AGN, or 
13 of 18 AGN, if Arp 220 templates are considered as 
starbursts. With either interpretation, the majority of 
sources are AGN. The comparable numbers for the H05 
sample were 30 of 31 AGN, or 27 of 31 AGN. The median 
24 /im flux density of the H05 sample is 1.1 mJy and of 
our present sample is 1.3 mJy, so both samples are very 
similar in this criterion. The current sample has median 
R > 24 mag (because 11 of 18 do not have measurable 
R, and the faintest measured is R = 24.5), whereas the 
H05 sample has median R > 25 mag (because the opti- 
cal survey limit in Bootes is fainter than in the FLS), 
so both samples are also similar in optical magnitude se- 
lection. The primary difference is the use of radio flux 
densities to define the current sample. Because the sam- 
ple was initially defined by the infrared limit, there was 
no constraint on the radio flux density as long as there 
was a radio detection. The radio sample reaches suffi- 
ciently faint to find sources with q as large as 1.0. The 
mean q for the 18 sources is 0.56. Only one source has 
a negative q; this is one of the four power law sources 
so is consistent with being an AGN. The median q value 
is a consequence of the fact that the great majority of 
sources in the FLS detected both at 24 /im and at 20 cm 
have approximately this value of q, a value co nsistent 
with that from starbursts l)Appleton et al. 120041) . There 
are numerous radio sources in the FLS with more nega- 
tive q, as expected from AGN, but these sources are not 
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found in large numbers in a sample requiring detection 
at 24 /im above flux densities of about one mJy. Because 
of the flux density limits used for our selection, our cri- 
teria lead, therefore, to a sample of objects which should 
be dominated by dusty, obscured starbursts according to 
our current understanding of the infrared and radio char- 
acteristics of such starbursts. As a result, we expected 
starburst spectral classifications to dominate our sample. 

In this context, it is significant that only a minority 
of sources show the starburst classification determined 
from the presence of PAH features. The mean q for the 
13 AGN (counting the power law sources), as classified 
in the IRS spectra, is 0.5; the mean for the 5 starbursts 
(counting the Arp 220 templates) is 0.6. Though ob- 
viously based on limited statistics, this result indicates 
that the q value is not a reliable diagnostic of whether a 
dusty source will show the PAH features expected from 
starbursts or the absorption expected from AGN. This 
is not a surprising result if we use Markarian 231 as the 
prototype dusty, absorbed AGN. At z = 2, it would have 
observed 24 /im flux density of 0.6 mJy, based on the rest- 
frame flux density of 1.5 Jy at 8 /im in lWeedman et al~1 
( 2005), so this absorbed AGN would be comparable to 
the flux densities of the infrared sources in our sample. 
More importantly, Markarian 231 would have a q value 
of 0.56, based on the ratio in the rest fra me of f„(8 /im) 
compa red to the f„(6 cm) of 0.41 Jy from lBecker et ahl 
(|1991|) : this value of q agrees very well with the median 
of the sources in our sample. 

4. SUMMARY AND CONCLUSIONS 

We have observed a new sample of 18 Spitzer 24 /im 
sources with the Spitzer IRS, chosen to be optically faint 
(R > 24 mag) and also to have radio detections at 20 
cm. The most definitive result of this study is confirma- 
tion of the conclusions of H05 that optically-faint Spitzer 
sources of f„(24/im) ~ 1 mJy are systematically at high 
redshift. Of the 18 sources in our sample, 12 (67 %) have 
confident redshifts z > 1.8. Within H05, at least 15 of 
31 (48 %) have z > 1.8. This is extraordinarily differ- 
ent from the redshift distribution expected from a flux 
limited 24 /im sample as predicted using current models 
or observations for sta rburst galaxies. For example, the 
redshift distribution bv lPerez- Gonzalez et al. 1 1)20051) for 
Spitzer sources to f„(24/im) ~ 0.1 mJy, based primar- 
ily on photometric redshifts for starburst templates, has 
< 10% of s o urces with z > 1.8. The earlier models of 
iDole et al. I l)2003fl predicted that for 24 /im samples to 
1.5 mJy, less than 1% would have z > 1.8. 

The redshift distribution of our optically-faint samples 
is also very different from the r esults for optically-br ight 
samples of similar infrared flux. Brown et al. I l)2005l) use 
a sample similar to ours in infrared flux limit (f„(24/im) 
> lmJy) but much brighter optically (R < 21.7). They 
have optically-determined redshifts of 255 sources in the 
Bootes field for which 61, or 24%, have z > 1.8, compared 
to the 67% which we find. Their selection is deliber- 
ately directed toward finding type 1 quasars by selecting 
objects with compact optical morphologies, and the re- 
sulting luminosity functions are similar to those derived 
from optically-selected quasar surveys to similar optical 
limits. It is clear, therefore, that an extreme infrared to 
optical flux ratio is a simple but powerful criterion for 
preferentially selecting dusty sources at high redshift. 



More challenging than determining the redshift is the 
question of classifying these dusty, high redshift sources 
as having their mid-infrared luminosities derived primar- 
ily from starbursts or AGN. This classification is essential 
for utilizing counts and redshifts of mid-infrared sources 
to determine the correct relative contribution of these 
two fundamental sources of luminosity in the early uni- 
verse. Our working hypothesis for this classification, 
based on analogies to local ultraluminous infrared galax- 
ies, is that sources showing strong absorption features 
or showing power-law continua without spectral features 
are classified as AGN, and those sources with conspic- 
uous PAH features are classified as starbursts. Using 
this initial hypothesis, we can examine some of the con- 
sequences of our results and make predictions regarding 
future results. Considering the results in Table 1, 16 
of the 18 sources show absorption features or power law 
spectra. For H05, the fraction is 30 of 31. In either sam- 
ple, this implies that over 90% of the sources are AGN. 
We can determine a lower limit to the AGN fraction by 
allowing that sources fit by the Arp 220 template are 
powered primarily by starbursts, because this template 
does include PAH emission even though it is dominated 
by absorption. In this case, 13 of the 18 sources in Table 
1 are AGN, and 27 of the 31 in H05. These diagnostics 
using IRS spectra indicate, therefore, that a large frac- 
tion of optically-faint sources in samples chosen at 24 /im 
are powered by AGN. Furthermore, these diagnostics in- 
dicate that at least 75% of those faint radio sources which 
would be attributed to star formation based on their in- 
frared to radio flux ratios are actually powered primarily 
by AGN. These results indicate that the AGN popula- 
tion at high redshift is higher than currently realized, 
in both radio and infrared samples, primarily because of 
a large fraction of optically obscured AGN which were 
previously unknown. 

There is a straightforward interpretation of these re- 
sults that explains comparisons to existing models of 
24 /im source counts and makes some testable predic- 
tions. The population of sources we have observed has 
f„(24/im) > 1 mJy and has extreme infrared to optical 
ratios such that R > 24 mag. If we conclude, as ex- 
plained above, that this population is primarily obscured 
AGN, it means that the population of high redshift, ob- 
scured AGN is much larger in number than the popu- 
lation of starbursts at mid-infrared luminosities corre- 
sponding to this flux density limit. It is not inconsistent, 
therefore, that the models of counts and redshift distri- 
butions that accomodate starbursts have not predicted 
this population of high luminosity, high redshift sources, 
because such models do not include the AGN popula- 
tion. This conclusion is also consistent with observations 
of the only known population of luminous, high redshift, 
dusty starbursts, which is the SCUBA or MAM BO sub- 
millimeter population l)Chapman et al. 1 120051) . These 
submillimete r sources are known to have f t/ (24/im) 
< 0-5 mJy l|Charmandaris et al. I l)2004t) iFraver et al. I 
(2004)). MAMBO observations of optically obscured 
sources selected to have f l/ (24/tm) > lmJy also indicate 
that these sources do not have the submillimeter fluxes 
expected for starbur st galaxies, but have relatively larger 
mid-infrared fluxes l)Lutz et al. 1120051?) . These various 
observations indicate that the luminous, dusty starburst 
population at z ~ 2 is not sufficiently luminous in the 
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mid-infrared to dominate the samples we selected for 
IRS observation. The important conclusion is that the 
high end of the infrared luminosity function for optically- 
obscured, dusty sources at z ~ 2 is dominated by AGN 
but that the relative fraction of starbursts increases for 
f„(24^m) < 0.5 mJy. This leads to a prediction about 
expected results from future observations of optically- 
faint sources selected only on the basis of f l/ (24^m): tar- 
gets observed at fainter 24 /im flux density levels should 
show a higher fraction of PAH-dominated spectra than 
in the samples we have presented. The one such source 
so far observed with f„(24/im) ~ 0.5 mJy, although cho- 
sen because of its submillimeter detection, does show 
PAH emission ijLutz et al. 112005^ . This conclusion also 
predicts that source characteristics as determined from 
overall spectral templates which fit a wide range of wave- 
lengths should indicate an increasing proportion of star- 
burst sources with decreasing f„(24/im), but populations 
with f„(24/xm) > 1 mJy should be dominated by AGN 
templates. 

If we predict that the nature of optically faint 24 /zm 
sources changes for f„(24/im) < 0.5 mJy (primarily star- 
bursts) compared to £,(24^m) > 1.0 mJy (primarily 
AGN), this also implies consequences regarding X-ray 
characteristics. In particular, comparison of X-ray and 
infrared detections should show a difference in the typ- 
ical infrared to X-ray flux ratio at these different in- 
frared flux limits. The mean infrared to X-ray flux ratio 
should be smaller for the brighter infrared sources than 
for the fainter infrared sources, because strong X-ray 
sources correspond to AGN. While compari sons of Chan- 
dra a nd Spitzer surveys have been made l)Rigbv et al. I 
I2004T) , these have been deep surveys in small areas of the 
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sky. As a result, there are insufficient sources in com- 
mon at the relevant flux levels to test this prediction 
meaningfully. For example, there are only 10 sources 
with X-ray detections an d f, y (24/nm) > l.QmJy i n thes e 
samples as displayed by lAlonso-Herrero et al. I f 2004) , 
and we do not know how many of those are optically 
faint or at high redshift. Our classification has defined 
sources with an absorbed infrared spectrum or power- law 
infrared spectrum as AGN. If this classification is correct, 
optically faint sources having f„(24/im) > lmJy chosen 
using X-ray selection as an additional criterion should 
always show power-law or absorbed spectra and never 
show PAH spectra, because the X-ray criterion is a firm 
AGN indicator. So far, no IRS results are available for 
sources chosen from combined Spitzer and Chandra sur- 
veys, but we have programs underway to obtain these. 
We can be optimistic, therefore, that as spectroscopic 
samples of Spitzer sources accumulate, we will be able to 
produce quantitative luminosity functions for both dusty 
starbursts and dusty AGN at the crucial epoch of z ~ 2. 
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Infrared Rcdshifts of Faint Radio Sources 



TABLE 1 

Properties of Sources 



Source 


Name 1 


f„ (24/im) 


R 


q b 


time c 


z(f) d 


z(t) c 






(mjy) 


(mag) 




(sec) 


1. 


SST24 J171440.57+584104.9 


1.30 




0.5 


1200 


1.98(abs)) 


2.00(231) 


2. 


SST24 J172210.27+584617.2 


2.28 


23.9 


1.0 


720 


1.16(abs) 


1.20(231) 


3. 


SST24 J172540.26+584953.4 


1.47 


24.2 


1.0 


720 


1.94(em) 


1.80(M82) 


4. 


SST24 J172249.03+585027.6 


1.12 


24.5 


0.7 


1440 






5. 


SST24 J171151.01+585104.1 


1.05 


24.4 


0.5 


1440 


1.94(abs) 


1.92(183) 


6. 


SST24 J171144.26+585225.8 


1.29 


23.9 


0.9 


720 


1.25(abs) 


1.22(231) 


7. 


SSI 24 J 171527.08+585802.2 


1.04 




-0.6 


1440 






8. 


SST24 J171205.09+591508.6 


1.10 




1.0 


1440 


2.31(abs) 


2.40(231) 


9. 


SST24 J171927.28+591536.5 


2.01 




0.4 


720 


2.25(em) 


2.20(M82) 


10. 


SST24 J170939.20+592728.3 


1.20 




0.7 


1440 


2.50(abs) 


2.52(220) 


11. 


SST24 J172708.53+593727.7 


1.11 




0.8 


1440 


1.96(abs) 


2.00(231) 


12. 


SST24 J171342. 79+593920.3 


1.89 




0.3 


720 






13. 


SST24 J172048.02+594320.6 


1.54 




0.3 


1200 


2.21(abs) 


2.29(183) 


14. 


SST24 J171343.94+595714.4 


1.03 


24.0 


0.4 


1440 


1.70(abs) 


1.81(183) 


15. 


SST24 J171057.45+600745.2 


2.05 


24.5 


0.7 


720 


2.38(abs) 


2.40(220) 


16. 


SST24 J172448.65+601439.1 


1.68 




0.4 


720 


2.22(abs) 


2.40(220) 


17. 


SST24 J172001.44+602544.9 


1.55 




0.2 


1200 


18. 


SST24 J171434.56+602828.6 


1.51 


24.2 


0.1 


1200 


2.13(abs) 


2.35(183) 



a SST24 source name derives from discovery with the MIPS 24/im images in the Spitzer FLS; coordinates listed are J2000 24/im positions 
with typical 3cr uncertainty of ± 1.2". Source 1 is extended in the 24/im image; the total flux density is about 3 mjy compared to the 
unresolved flux density of 1.3 mjy listed in the Table. 

b q = log[f„(24 (um)/f„(20 cm)] 

Tntegration time for each order of Long Low spectrum; integration time in Short Low order 1 is 240 s in all cases. 

d rcdshift z(f) is the redshift that would be assigned by locating apparent features in the spectrum and identifying these with PAH emission 
features (em) at 6.2ttm, 7.7/^m and 11.2/im (rest frame) or the continuum excess centered at 8fim arising just shortward of the 9.7/im 
silicate absorption feature (abs). 

°redshift z(t) is the redshift determined from a formal chi-squared template fit as illustrated in Figures 1, 2 and 3; templates are identified 
in parenthesis as pure PAH emission (M82); combination of PAH emission and silicate absorption as seen in Arp 220 (220); pure absorption 
as seen in Markarian 231 (231), and very deep absorption as seen in IRAS F00183-7111 (183). We note that source 13 has been listed by 
Martincz-Sansigre et al. (2005) as having a weak Ly a line at z = 2.02. 

t a is the power law index for the featureless continuum in sources without measurable redshift. 
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Fig. 1. — Observed spectra smoothed to approximate resolution of individual IRS orders (histogram), best template fit (solid line), 
and derived rcdshift for all sources in Table 1 with unambiguous fit for redshift or with no fit; panels without a template identification or 
rcdshift are sources for which there are insufficient spectral features to derive a redshift; panels arc labeled by number in Table 1. 
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Fig. 2. — Observed spectra smoothed to approximate resolution of individual IRS orders (histogram) and alternative template fits (solid 
lines) from either absorption templates or PAH emission templates and the resulting rcdshifts for sources for which either silicate absorption 
or PAH emission templates could provide a possible fit; panels are labeled by number in Table 1; first panel shown for each source is the 
fit selected and included in Table 1. 
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Fig. 3. — Observed spectra smoothed to approximate resolution of individual IRS orders (histogram) and alternative template fits (solid 
lines) from cither absorption or PAH emission and the resulting redshifts for sources for which cither silicate absorption or PAH emission 
templates could provide a possible fit; panels are labeled by number in Table 1; first panel shown for each source is the fit selected and 
included in Table 1 except for source 17, for which no redshift is assigned. 



